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The value of the Prandtl number P exerts a strong influence on convection driven dy- 
namos in rotating spherical shells filled with electrically conducting fluids. Low Prandtl 
numbers promote dynamo action through the shear provided by differential rotation, 
while the generation of magnetic fields is more difficult to sustain in high Prandtl num- 
ber fluids where higher values of the magnetic Prandtl number P m are required. The 
magnetostrophic approximation often used in dynamo theory appears to be valid only 
for relatively high values of P and P m . Dynamos with a minimum value of P m seem 
to be most readily realizable in the presence of convection columns at moderately low 
values of P. The structure of the magnetic field varies strongly with P in that dynamos 
with a strong axial dipole field are found for high values of P while the energy of this 
component is exceeded by that of the axisymmetric toroidal field and by that of the 
non-axisymmetric components at low values of P. Some conclusions are discussed in re- 
lation to the problem of the generation of planetary magnetic fields by motions in their 
electrically conducting liquid cores. 



1. Introduction 

The problem of the generation of magnetic fields by motions of an electrically con- 
ducting fluid in rotating spherical shells is one of the fundamental problems of planetary 
and astrophysical sciences. The increasing availability in recent years of computer ca- 
pacity has permitted large scale numerical simulations of this process, but the creation 
of realistic models for planetary and stellar dynamos has been hampered by a lack of 
knowledge about appropriate external parameters. Since molecular values of material 
properties are usually not attainable in computer simulations because of the limited nu- 
merical resolution, eddy diffusivities must be invoked for comparisons with observations. 
Eddy diffusivities represent the effects of the unresolved scales of the turbulent velocity 
field and it is often assumed for this reason that the eddy diffusivities for velocities, 
temperature and magnetic fields are identical. The effects of turbulence on the diffusion 
of vector and scalar quantities differ, however, and large differences in the corresponding 
molecular diffusivities are likely to be reflected in the effective diffusivities caused by the 
fluctuating velocity field of the unresolved scales. In the case of the Earth's core, for 
example, the magnetic diffusivity is assumed to have a value of the order of 2 m 2 / sec 
(Braginsky and Roberts, 1995) which exceeds the kinematic viscosity by a factor of the 
order 10 6 . Numerical simulations are thus capable of resolving magnetic fields, but are 
far from resolving velocity structures. However, since the largest unresolved scales, v, I, 
of velocity and length yield values for an eddy viscosity v e ss v ■ I of the order 1 m 2 / sec 
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Figure 1. Geometrical configuration of the problem. A part of the outer spherical surface is 
removed to expose the interior of the shell to which the conducting fluid is confined. 

or less, it is reasonable to assume magnetic Prandtl numbers less than unity. Since the 
concept of eddy diffusivities is simplistic there have not been many theoretical deriva- 
tions for ratios of turbulent diffusivities which enter the dimensionless equations for the 
numerical simulations. Some theoretical considerations in the astrophysical context can 
be found in the paper by Eschrich and Riidiger (1983). Values for an effective Prandtl 
number can eventually be derived from experiments on turbulent convection (see, for 
example, Ahlers and Xu, 2001). Additional complications arise in planetary and astro- 
physical applications through anisotropies introduced by the effects of rotation and the 
presence of large scale magnetic fields. While anisotropic eddy diffusivity tensors will 
be useful for more realistic models, for the purpose of the present paper we prefer the 
simplicity of scalar diffusivities. 

After introducing the basic equations and the method of their numerical solution in 
section 2 we shall consider the most important properties of convection without magnetic 
field as a function of the Prandtl number in section 3. Since a number of earlier papers on 
this topic have been published this section can be kept relatively short. In section 4 the 
onset of convection driven dynamos in fluids with different Prandtl numbers is described. 
Energy aspects are considered in section 5 and the validity of the magnetostrophic ap- 
proximation is discussed in section 6. The influences of various boundary conditions are 
considered in section 7 and a concluding discussion is given in the final section 8. 



2. Mathematical formulation of the problem and methods of solution 

We consider a rotating spherical fluid shell as shown in figure [TJ We assume that a 
static state exists with the temperature distribution Tg = T — (3d 2 r 2 /2. Here rd is the 
length of the position vector with respect to the center of the sphere. The gravity field is 
given by g = —djr. In addition to d, the time d 2 /v, the temperature v 2 /'yad 4 and the 
magnetic flux density v{\iq) x I 2 jd are used as scales for the dimensionless description of 
the problem where v denotes the kinematic viscosity of the fluid, k its thermal diffusivity, 
g its density and \i is its magnetic permeability. The equations of motion for the velocity 



Prandtl number dependence of convective dynamos 367 

vector u, the heat equation for the deviation O from the static temperature distribution, 
and the equation of induction for the magnetic flux density B arc thus given by 

d t u + u ■ Vu + rfc x u = -Vtt + 9r + V 2 w + B ■ VB, (2.1a) 

V-it = 0, (2.16) 

P(d t e + u ■ V6) = Rr ■ u + V 2 6, (2.1c) 

V-B = 0, (2.1rf) 

V 2 B = P m {d t B + u-VB - B - Vu), (2.1 e) 

where d t denotes the partial derivative with respect to time t and where all terms in 
the equation of motion that can be written as gradients have been combined into V7r. 
The Boussinesq approximation has been assumed in that the density g is regarded as 
constant except in the gravity term where its temperature dependence given by a = 
— (dg/ dT)/g —const is taken into account. The Rayleigh number R, the Coriolis number 
r, the Prandtl number P and the magnetic Prandtl number P m are defined by 

R=^—, r= , P=~, P m = ~, (2.2) 

VK V K A 

where A is the magnetic diffusivity. Because the velocity field u as well as the magnetic 
flux density B are solenoidal vector fields, the general representation in terms of poloidal 
and toroidal components can be used 

u = V x (V« x r) + Ww x r , (2.3a) 
B = V x (Vft x r) + V.g x r . (2.36) 

By multiplying the (curl) 2 and the curl of equation (|2.1a|) by r we obtain two equations 
for v and w 

[(V 2 - d t )C 2 + rd v ]\7 2 v + tQw - £ 2 6 = -r ■ V x [V x (u ■ Vu - B ■ VB)], (2.4a) 

[(V 2 - d t )C 2 + Td v ]w - tQv = r -V x (u-Vu- B ■ VB), (2.46) 

where d v denotes the partial derivative with respect to the angle ip of a spherical system 
of coordinates r, 9, tp and where the operators C 2 and Q are defined by 

£2 = -r 2 V 2 + d r (r 2 d r ), 

Q = rcos6V 2 - (£ 2 + rd r )(cos6d r - r" 1 smOdg). 

The heat equation for the dimensionless deviation from the static temperature distri- 
bution can be written in the form 

V 2 e + R£ 2 v = P{d t +u- V)9, (2.5) 

and the equations for h and g are obtained through the multiplication of equation (|2.1e| 
and of its curl by r 

V 2 C 2 h = P m [d t C 2 h -fVx(uxB)], (2.6a) 
V 2 C 2 g = P m [d t C 2 g - r ■ V x (V x (u x B))}. (2.66) 

We assume stress-free boundaries with fixed temperatures and use the value 0.4 for the 
radius ratio rj = r,/r , 

v = d 2 r v = d r (w/r) = 9 = at r = n = 2/3 and r = r Q = 5/3. (2.7) 

For the magnetic field electrically insulating boundaries are assumed such that the 
poloidal function h must be matched to the function which describes the poten- 



368 



R. Simitev and F. H. Busse 




Figure 2. Equatorial streamlines rd v v = const, in the case P = 0.5, r = 10 4 and R = 1.8 X 10 5 
3.2 x 10 5 , 4 x 10 5 , 9 x 10 5 (from left to right). The two halves of last plot show convection in 
the minimum and in the maximum of a relaxation cycle. 



tial fields outside the fluid shell 

g = h-h {e) = d r {h-h {e) ) = at r = r t = 2/3 and r = r Q = 5/3. (2.8) 

But computations for the case of an inner boundary with no-slip conditions and an 
electrical conductivity equal to that of the fluid have also been done. The numerical 
integration of equations (|2.4|) . (|2.5|) and (|2.6[) together with boundary conditions (|2.7p 
and (|2.8p proceeds with the pseudo-spectral method as described by Tilgner and Busse 
(1997) which is based on an expansion of all dependent variables in spherical harmonics 
for the 9, ^-dependences, i.e. 

v = Vi m (r, t)Pr {cos 8) exp{imp} (2.9) 

Lrn 

and analogous expressions for the other variables, w, 0, h and g. P[™ denotes the asso- 
ciated Legendre functions. For the r-dependence expansions in Chebychev polynomials 
are used. For further details see also Busse et al. (1998) or Grote et al. (1999). For the 
computations to be reported in the following a minimum of 33 collocation points in the 
radial direction and spherical harmonics up to the order 64 have been used. But in many 
cases the resolution has been increased to 49 collocation points and spherical harmonics 
up to the order 96 or 128. 

3. Convection in rotating spherical shells 

Three different types of convection can be distinguished in rotating spherical shells. 
Predominantly convection occurs in the form of rolls aligned with the axis of rotation 
which exhibit properties of thermal Rossby waves in that they are drifting in the prograde 
azimuthal direction. They are confined to the region outside the virtual surface of the 
tangent cylinder which touches the inner boundary at its equator. The dynamics of these 
convection columns, as they are sometimes called, is intimately connected for Prandtl 
numbers of the order unity or less with the differential rotation that is generated by 
their Reynolds stresses. Above their onset as a m-periodic pattern in the azimuthal 
direction the convection columns experience with increasing Rayleigh number transitions 
to amplitude and shape vacillations before they become spatio-temporally chaotic in 
the dimensions perpendicular to the axis while retaining their nearly perfect alignment 
with the rotation vector. In this regime of beginning turbulence coherent processes such 
as localized convection and relaxation oscillations are realized through the interaction 
of the differential rotation and the convection columns. A graphical display of these 
stages of convection is shown in figure [2] in terms of the streamlines of the convection 



Prandtl number dependence of convective dynamos 369 




Figure 3. Equatorial streamlines rd v v = const, in the case P = 20, r = 5 x 10 3 and 
R = 1.6 X 10 5 , 1.75 x 10 5 , 2.5 x 10 5 , 1.5 x 10 6 (from left to right). 




Figure 4. Equatorial streamlines rd v v = const, in the case P = 0.025, t = 10 5 and R = 4 x 10 5 , 
6 x 10 5 , 8 x 10 5 , 10 6 (from left to right). 



columns intersected by the equatorial plane. For a more detailed description we refer 
to Grote and Busse (2001) or the review of Busse (2002a). As the Reynolds stresses 
decrease with increasing P these coherent processes disappear. Figure [3] indicates that 
the spiralling nature of the convection columns also diminishes with increasing Rayleigh 
number. A small non-geostrophic differential rotation persist driven as a "thermal wind" 
by temperature gradients caused by the lateral inhomogeneity of the convective heat 
transport. 

A second distinct form of convection are the equatorially attached cells which represent 
modified inertial modes and become the preferred form of convection at sufficiently low 
values of P. They were first found by Zhang and Busse (1987) and analytical descrip- 
tions in terms of perturbed inertial oscillations have been given by Zhang (1994, 1995) 
and Busse and Simitev (2004a). The equatorially attached cells do not develop strong 
Reynolds stresses and they are thus less subject to the disruptive effects of the shear of 
a differential rotation. But the continuity of the convective heat transport requires that 
the equatorially attached convection occurs in conjunction with the columnar convection 
closer to the inner boundary of the shell. This effect is visible in the form of secondary 
extrema of the streamlines far from the boundary in figure [?] where the evolution of the 
convection flows in the chaotic regime is illustrated. Note that the strong attachment to 
the outer equatorial boundary persists even at the highest value of R used in this figure. 
For further details see Simitev and Busse (2003). 

A third form of convection is realized in the polar regions of the shell which comprise 
the two fluid domains inside the tangent cylinder. Since gravity and rotation vectors are 
nearly parallel in these regions (unless large values of rj are used) convection resembles 
the kind realized in a horizontal layer heated from below and rotating about a vertical 
axis. A tendency towards an alignment of convection rolls with the North-South direction 
(Busse and Cuong, 1977) can be noticed, but this property is superseded by instabilities 
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Figure 5. Polar convection in the cases, P = 0.025, r = 5 x 10 4 , R = 2 x 10 5 ; P = 0.1, 
r = 3 x 10 4 , R = 7.5 x 10 5 ; P = 1, r = 10 4 , P = 1.4 x 10 6 and P = 20, r = 5 x 10 3 , 
P = 10 6 (from left to right). The plots in the upper row show lines of constant u T at the surface 
r — n + 0.5 as seen from the north pole. Isolines corresponding to higher amplitude have not 
been plotted in order to emphasize structuresjvith lower amplitude. The plots in the lower row 
exhibit lines of constant u v in the left half, of O in the upper right quarter and of the streamlines 
rsmOdgv in the lower right quarter, all in the meridional plane. 



of the Kiippers-Lortz type and by interactions with turbulent convection outside the 
tangent cylinder. The onset of convection in the polar regions generally occurs at Rayleigh 
numbers considerably above the critical values R c for onset of convection outside the 
tangent cylinder. Except for the case of very low Prandtl numbers the differential rotation 
in the polar regions is usually oriented in the direction opposite to that of rotation 
and thus tends to facilitate polar convection by reducing the rotational constraint. The 
possibility exists, however, that at sufficiently low values of P and high values of r finite 
amplitude convection in the polar regions may precede the onset of convection in other 
regions. The patterns of polar convection and other properties discussed in this section 
are illuminated by the display of convection at different values of P in figure [5] 

Quantitative aspects of convection are described by figure [6] where the averages over 
space and time of the kinetic energy densities of the various components of the convection 
flow are shown as a function of R. The energy densities are defined by 

1(| Vx (Vvxr) | 2 ), E t = 1(| Viixr | 2 ), (3.1a) 

~(| Vx (Vixr) | 2 ), Et = 1(| Vw x r | 2 ), (3.16) 

where the angular brackets indicate the average over the fluid shell and v refers to the 
azimuthally averaged component of v, while v is defined by v = v — v. The Nusselt 
number at the inner spherical boundary Nui is also shown in figure [SI It is defined by 

(3.2) 

where the double bar indicates the average over the spherical surface. The rapid growth 
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{R-R c )/R c (R~R C )/R C 

Figure 6. Time- averaged energy densities E t (squares), E p (plus-signs), E t (crosses) andNusselt 
number Nm (filled triangles, right ordinates) as functions of R/R c for (a) P — 0.025, r = 10 5 , 
(b) P = 0.1, t = 3 x 10 4 , (c) P = 1, t = 10 4 , and (d) P = 20, r = 5 x 10 3 . The values of R c 
are 283000, 235000, 190000 and 146000 in the cases (a), (b), (c) and (d) respectively. The mean 
poloidal energy component E p has been omitted since it is by factor 10 3 -10 2 smaller than the 
other components. 



with R of E t corresponding to the energy of differential rotation is remarkable in the cases 
(a), (b), and (c) of figure^ Only for higher values of P as in case (d) does E t never exceed 
the energies of the fluctuating components of motion. Another remarkable feature of the 
plots is the rapid rise of the Nusselt number Nm after the onset of amplitude vacillations 
and the transition to a chaotic time dependence. The steadily drifting convection columns 
or the equatorially-attached cells do not transport heat very well because of the mismatch 
between the geometry of the boundaries and the flow structure. Because this mismatch 
does not occur in the case of polar convection, the heat transport of the latter may 
easily exceed that of convection outside the tangent cylinder at higher values of R. The 
dependence of the kinetic energy densities on R in case (a) of figure [S] are somewhat 
unusual in that after a rapid growth near onset a regime of saturation follows and a 
growth of the kinetic energies with R begins again only when R has reached nearly 
twice its critical value. The behavior is caused by the competition of several modes with 
different azimuthal wavenumbers m as is described in more detail by Simitev and Busse 
(2003) . The ways in which properties of convection are reflected in the generated magnetic 
fields will be discussed in the following sections. 
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Figure 7. Magnetic Reynolds numbers Rm for the onset of dynamo action as a function of P m 
in the cases P = 0.01, r = 10 5 (stars), P = 0.025, r = 10 s (crosses), P = 0.1, r = 10 s (circles), 
P = 1, r = 3 x 10 4 (triangles up), P = 1, r = 10 4 (squares) and P = 5, r = 5 x 10 3 (diamonds) 
and P — 10, r = 5 x 10 3 (triangles down). The empty symbols in the cases with P ^ 0.1 are 
based on decaying dynamos, the full symbols are based on the lowest non-decaying solutions. 



4. The onset of dynamo action for different Prandtl numbers 

The onset of convection driven dynamos in the case P = 1 has been explored in nu- 
merous papers in the past. We refer to the work of Kageyama and Sato (1997), Kida and 
Kitauchi (1998), Ishihara and Kida (2002), Busse et al. (1998), Christensen et al. (1999), 
Olson et al. (1999), Katayama et al. (1999), Takahashi et al. (2003), Grote et al. (2000), 
Grote and Busse (2001). In particular the results of the latter two papers can be com- 
pared well with the results presented in the following since only the Prandtl number differs 
among the parameters of the computations. Since one of the goals of the numerical sim- 
ulations is the search for dynamos with low values of the magnetic Prandtl number P m , 
the amplitude of convection must be sufficiently high such that the magnetic Reynolds 
number defined by Rm — (2E) 1 ' 2 P m reaches a value of the order 10 2 . Here E refers to 
the average density of the kinetic energy of convection, i.e. E = E p + E t + E p + E t . 
In figure [7] values of Rm for cases of sustained dynamo action and for cases of decaying 
magnetic fields have been indicated for several values of P. No special effort has been 
made to determine the minimum value Rm m i n of Rm for dynamo action. Since in all 
cases convection is chaotic and since dynamos are typically subcritical, i.e. the boundary 
between sustenance of a dynamo and its decay occurs at finite amplitudes of the magnetic 
field, Rm m i n is not a well defined quantity. Moreover, the value Rm as defined above 
will decrease after the onset of dynamo action since typically the energy of differential 
rotation is strongly reduced by the action of the Lorentz force. 

It is evident from figure [7| that the goal of low values P m is most readily reached at 
low Prandtl numbers. At high values of P the minimum value of the magnetic Reynolds 
numbers for dynamos decreases slightly, but P m must increase in proportion to P in order 
to achieve a finite ratio P/ P m in the limit of infinite Prandtl number as shown by Zhang 
and Busse (1990). The tendency of low Prandtl numbers to lower the minimum values of 
P rn required for dynamos is limited, however, by the transition to inertial convection for 
P ^ (t/t )~ 2 where t is a quantity of the order 10 2 (Ardes et al., 1997). It thus has not 
been possible to obtain dynamos with values of P m below 0.1 for Prandtl numbers less 
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Figure 8. Convection driven dynamos as a function of R, P and P m for r = 5 x 10 3 . The symbols 
indicate chaotic dipolar (a), hemispherical (A), quadrupolar mixed (■) and decaying 

dynamos (©). 

than 0.05. Inertial convection is evidently less conducive to dynamo action than columnar 
convection and higher values of R m are required for dynamos at P = 0.025 and P = 0.01 
than at higher Prandtl numbers according to figure [7J The most promising way towards 
low P m dynamos thus appears to be an increasing r with an accompanying reduction of 
P. Of course, the numerical resolution will have to be increased at the same time. 

A particular feature that is evident from figure [7| (and also figure [9|) is that high values 
of R may prevent dynamo action. In the case P — 0.1, t — 10 5 it is found that only a 
finite interval of convection energy permits dynamo action since high values of Rm lead 
to the expulsion of magnetic flux from the convection eddies which is detrimental for the 
generation of magnetic fields. The same effect can also be noticed for higher values of 
the Prandtl number as demonstrated in the case of P m = P = 10 of figured! 

A more complete overview of dynamos at values of P and P m above unity is shown 
in figure [8] where the structure of the magnetic field is also indicated. As long as the 
properties of convection are nearly symmetric with respect to the equatorial plane three 
types of dynamos can be distinguished: Quadrupolar dynamos are characterized by a 
magnetic field that exhibits the same symmetry with respect to the equatorial plane as 
the convection velocity field, while the opposite symmetry characterizes dipolar dynamos. 
Hemispherical dynamos correspond to the superposition of a quadrupolar and a dipolar 
magnetic field of nearly the same form, but with opposite signs in the northern and 
the southern hemispheres such that the magnetic field nearly vanishes in one of the 
hemispheres (Grote and Busse, 2000). Once the symmetry of convection with respect 
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Figure 9. Convection driven dynamos as a function of the Rayleigh number R and the magnetic 
Prandtl number P m for P — 0.1, r = 10 5 (empty symbols) and r = 3 x 10 4 (filled symbols). The 
symbols indicate chaotic dipolar (squares) , hemispherical (triangles) , mixed (stars) and decaying 
dynamos (circles). 



to the equatorial plane is impaired significantly, - as happens, for example, at higher 
Rayleigh numbers through the onset of polar convection -, the structures of dynamos 
can no longer be clearly distinguished and we speak of "mixed" dynamos. These dynamos 
continue to exhibit an alignment with the axis of rotation. Dynamos with a dominant 
equatorial dipole, for instance, have not been found. 

For a fixed value of P quadrupolar dynamos are obtained for sufficiently low values 
of P m which give way to hemispherical and dipolar dynamos with increasing P m as has 
also been observed in the case P = 1 (Grote et al., 2000). But since dynamos for low 
values of P m disappear with increasing P quadrupolar and hemispherical dynamos can 
no longer be obtained as P exceeds a value of the order 5 in the case of r = 5 x 10 3 . 
Quadrupolar dynamos thus appear to be restricted to Prandtl numbers of the order unity 
and moderate values of r. When r is increased from 10 4 to 3 x 10 4 quadrupolar dynamos 
do no longer seem to be accessible according to the results of Grote et al. (2001, see also 
Busse, 2002a). The same effect occurs as P is lowered as indicated in figure M for the case 
P = 0.1, r = 10 5 . Here the lowest value of P m corresponds to a hemispherical dynamo. 

It should be noted that the critical value R c of the Rayleigh number increases nearly 
in proportion to rP for P < 1 and nominally higher values of r can thus be reached for 
lower values of P. The degree of chaos increases, however, with decreasing P because of 
the short thermal time scale. Nevertheless coherent structures are also found such as the 
relaxation oscillations apparent in the time record shown in figure 1101 Here the energy 
densities on the various components of the magnetic field are shown which are defined 
in analogy to the kinetic energy densities (|3.ip , 

Mp^dVxfVfexr)! 2 ), ~M t = i(| V<? x r | 2 ), (4.1a) 

M P = i<| V x (Vh x r) | 2 ), M t = i(| V<? x r | 2 ). (4.16) 

In contrast to the relaxation oscillations of non-magnetic convection mentioned in con- 
nection with figureOthe oscillations of figure [TD] have a much shorter period and originate 
from the oscillations of the hemispherical magnetic field as is evident from the comparable 
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Figure 10. Time series of energy densities of a hemispherical dynamo in the case P = 0.1, 
r = 10 5 , R = 6x 10 6 , P m = 0.11. The upper and middle panels show energy densities of dipolar 
and quadrupolar components of the magnetic field, while the lower panel shows energy densities 
of the velocity field. The mean toroidal components are represented by solid lines, the fluctuating 
toroidal - by dotted lines, the mean poloidal - by dot-dashed lines and the fluctuating poloidal 
by dashed lines. 

periods. Since the strength of the magnetic field varies more with the phase of the cycle 
at low values of P than at higher ones the action of the Lorentz force on the differential 
rotation also exhibits significant variations. This in turn affects the amplitude of con- 
vection and its dynamo action and leads to the cyclic feedback process evident in figure 
[TU1 We note that the computational results displayed in figure [TU] were obtained with a 
resolution of 33 collocation points in the radial direction and with spherical harmonics 
of the order 96. Results obtained with 41 collocation points and spherical harmonics up 
to the order 108 are rather similar and the time averaged energies differ insignificantly 
Finally, we like to point out some typical differences between high and low Prandtl 
number dipolar dynamos which are exhibited by the plots of figure 1111 At low values 
of P the quadrupolar component is usually quite noticeable leading to a shift of the 
magnetic equator away from the geometric one. Often this shift alternates in time in 
a nearly periodic fashion. In high Prandtl number dynamos the mean magnetic field is 
nearly steady and the quadrupolar component is negligible if the Rayleigh number is not 
too high. This is evident in the lower plots of figure QT] even though the onset of polar 
convection with its attendant asymmetry with respect to the equatorial plane can already 
be noticed. A characteristic feature of high P dynamos are the strong zonal magnetic 
flux tubes in the polar regions which are mainly caused by the thermal wind shear an 
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Figure 11. Non-oscillatory chaotic dipolar dynamos in the cases P = 0.1, r = 10 5 , R = 2x 10 6 , 
P m = 3 (upper row) and P = 200, r = 5 x 10 3 , R = 10 6 , P m = 80 (lower row). The left column 
shows meridional isolines of B v (left half) and of rsmOdeh (right half). The middle column 
shows lines B r = const, at r — r<j. The right column shows lines u r — const, at r — + 0.5. 

example of which can be seen in the case P = 20 of figure 5. These flux tubes are also 
very evident in the dynamo of Glatzmaier and Roberts (1995). While the differential 
rotation changes a lot with increasing Prandtl number, the structure of the convection 
columns hardly varies as is also evident from figure [TT] The rate of drift in the azimuthal 
direction decreases, of course, with increasing P. 



5. Energetics of dynamos 

It is of interest to compare the various interactions between velocity and magnetic 
field components which sustain dynamo action against Ohmic dissipation. Zhang and 
Busse (1989) have listed a total of 31 interaction integrals which contribute to the energy 
balances obtained from equations (|2.6p for the mean and fluctuating components of the 
poloidal and toroidal magnetic fields. Fortunately only a few of the 31 terms contribute 
significantly and we list here only the most important ones, 

p 3 = (vhh), (5.1a) 

i 3 = (vgg), (5.16) 

p 3 ee (vhh), (5.1c) 

t 3 ee (whg), (5.1d) 

where the first two letters inside the brackets indicate which of the interactions between 
velocity and magnetic field components on the right hand sides of equations (|2.6p coun- 
teract the Ohmic dissipation of the magnetic field component indicated by the last letter 
inside the brackets. In the case of chaotic dynamos these integrals tend to fluctuate 
wildly and may change their signs. The particular ones listed above have been plotted as 
a function of time in figure [T^] together with the corresponding average Ohmic dissipa- 
tion density for three representative cases. It is remarkable that the differential rotation 
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Figure 12. Lorentz terms (the thin solid, dotted, and dashed lines indicate the first, second, 
and third terms, respectively, in each of the lines of (5.1)) and Ohmic dissipations (thick dotted 
shaded lines) in the cases P = 0.1, r = 10 5 , R = 6 x 10 6 , P m = 0.15 (left column), P = 0.1, 
t = 10 5 , R = 3 x 10 6 , P m = 3 (middle column), P = 200, r = 5 x 10 3 , R = 10 6 , P m = 80 (right 
column). 

contributes most of the sustenance of the mean toroidal field not only in the case of 
low values of P and P m , but also in the high Prandtl number case. The mean thermal 
wind is obviously sufficient to generate the mean toroidal field through the distortion 
and stretching of the mean poloidal field. This w-effect is thus usually found to dominate 
the dynamo process in the case of the zonal field. Only in particularly chaotic dynamos 
characterized by a rather weak mean poloidal field is the w-effect not effective as in the 
case P = 0.1, P m = 3 of figure [T2J ft is remarkable that the first and second terms of 
definitions (5.1a) are strongly anticorrelated in the low Prandtl number cases of figure 
1121 Because of their near cancellations the axisymmetric poloidal field is rather weak in 
low P dynamos. 

The a-effect in which the fluctuating magnetic field is generated through the interaction 
between fluctuating velocity field and mean magnetic field plays a far less important role 
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Figure 13. Kinetic E x and magnetic M x energy densities (left) and viscous V x and Ohmic O x 
dissipations (right) as functions of P m for convection driven dynamos for r = 5 x 10 3 , R = 10 6 
and Prandtl number as indicated in the boxes. The highest values of the Elsasser number for 
the cases P = 1, 3, 5, 15 are A = 3.02, 0.31, 0.37 and 0.49, respectively. The components X p , 
Xt, X p , Xt (where X = E, M, V, O) are represented by circles, squares, plus-signs and crosses, 
respectively. Kinetic energy densities and viscous dissipations are shown with light symbols, 
magnetic energy densities and Ohmic dissipations are shown with heavy symbols. 



than the w-effect. The fluctuating components of the magnetic field are usually generated 
through interactions of fluctuating parts of velocity and magnetic fields except in the case 
of high Prandtl numbers where the mean poloidal field typically dominates. 

One of the least understood questions of convection driven dynamos is the equilibration 
of magnetic energy. While in the astrophysical context an equipartition between kinetic 
and magnetic energy is often favored, an Elsasser number A of the order unity is regarded 
by geophysicists as a good estimate of the magnetic energy generated by dynamos in 
planetary cores. In the following we shall use the definition 

A=^ (5.2) 

T 

for the Elsasser number where M denotes the magnetic energy density averaged over the 
fluid shell and in time, i.e. M = M p +JI t + M p + M t . Already Chandrasekhar (1961) 
established that for A = 1 a minimum is found of the critical Rayleigh number for onset 
of convection in a horizontal layer heated from below, penetrated by a homogeneous 
vertical magnetic field and rotating about a vertical axis. 

In figures [13] and [HI averaged kinetic and magnetic energy densities and viscous and 
ohmic dissipation densities have been plotted for typical parameter sets of the problem. 
The Elsasser number has been computed for the largest magnetic energy M in each of 
the boxes and is given in the figure captions. It is evident that A rarely exceeds unity 
and when it does a physical reason of a high value of A is not obvious. In particular, in 
the cases of figure [13] where only P m is varied it is hard to understand why the Elsasser 
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Figure 14. Kinetic E x and magnetic M x energy densities (upper row) and viscous V x and 
Ohmic O x dissipations (lower row) as functions of R for convection driven dynamos in the cases 
(a) P = 0.1, r = 10 5 , P m = 0.2, A = 0.38; (b) P = 0.1, r = 10 5 , P m = 3, A = 3.22; (c) 
P = 0.1, r = 3 x 10 4 , P m = 1, A = 1.87; (d) P = 1, r = 3 x 10 4 , P m = 2, A = 0.26; (e) P = 1, 
r = 5 x 10 3 , P m = 1, A = 0.93; (f) P = 5, r = 5 x_10 3 J _P m = 10, A = 0.68 and (g) P = 10, 
r = 5 x 10 3 , P m = 10, A = 0.12. The components X p , X t , X p , X t (where X = E, M, V, O) 
are represented by circles, squares, plus-signs and crosses, respectively. Kinetic energy densities 
and viscous dissipations are shown with light symbols, magnetic energy densities and Ohmic 
dissipations are shown with heavy symbols. In case (b) results for a convection driven dynamo 
as well as for convection without magnetic field have been plotted in the case R — 2 x 10 6 . The 
values of A represent the highest values for each of the boxes. 
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Figure 15. Convection driven dynamo for P = 0.01, r = 10 5 , R = 6 x 10 5 and P m = 0.5. The 
first plot shows lines of constant u v in the left half and streamlines, rsmOdev = const, in the 
right half. The second plot shows streamlines in the equatorial plane, rd v v — const.. The left 
half of the third plot shows lines of constant zonal magnetic field and the right half shows 
meridional field lines, r sin 9dgh — const.. The last plot shows lines of constant B r at the surface 
r — r a + 0.5. 

number increases more strongly than P m itself. Another way of looking at the problem is 
to determine the highest value of A as a function of the Rayleigh number. A systematic 
study of the variation of magnetic energies with increasing R has been undertaken by 
Grote and Busse (2001) in the case P — P m = 1, t = 5 x 10 3 . This study has been 
extended since that time in order to obtain reliable time averages, see figure 11 of Busse 
and Simitev (2004b). Here the case of maximal M is attained for R = 1.4 x 10 6 in part 
(e) of figure [T4l and corresponds to a value of about unity for A. But this result may be 
accidental as is suggested by the magnetic fields in the cases of dynamos for low values 
of P and P m . Here again we find in the case (a) of figure [T4l that magnetic energy has 
reached a value close to its maximum as a function of i?, but the corresponding value 
0.38 of A is significantly lower. A much higher value of A is obtained, on the other hand, 
when P m = 0.2 is replaced by P m = 3 even though M may still increase with a further 
increase in R. We thus conclude that the criterion A w 1 for the equilibration of the 
magnetic energy is useful as a rough estimate, but the dependence of the equilibrium 
magnetic energy on P m is actually weaker than suggested by the definition (|5.2|l of A. 

There are a number of other features exhibited by the diagrams of figures Q2] and [HJ 
While the energy densities of the fluctuating components of the magnetic field usually 
exceed those of the mean components this situation is reversed at sufficiently high Prandtl 
numbers for the poloidal part of the field as is evident from the last columns of figure 
[TBI as well as of figure [Tj] This changeover occurs rather suddenly at about P = 8 for 
t = 5 x 10 3 as can be noticed in the comparison of cases (f) and (g) of figure [T4l A closer 
inspection shows that the replacement of the geostrophic differential rotation by the 
thermal wind and the accompanying growth of the polar zonal flux tubes is responsible 
for this change in the dynamo process. It should be noticed that this changeover is not 
evident in the plots of the dissipation densities since the latter are always much lower for 
the mean components than for the fluctuating components. 

Another property that changes with increasing Prandtl number is the fact that for 
P of the order unity or less the energy of the mean toroidal field exceeds that of the 
mean poloidal field except for some high Rayleigh number cases where polar convection 
begins to dominate. For dipolar dynamos obtained for P = 0.01, r = 10 5 with Rayleigh 
numbers of the order 10 6 and 0.2 ^ P m ^ 0.5 the energy of the mean toroidal field 
exceeds that of the mean poloidal field by more than a factor 10 2 . For example, in the 
case R = 10 6 ,P m = 0.5 the time averaged magnetic energy densities M p = 146 and 
M t — 2.26 x 10 4 are found. The structure of velocity and magnetic fields in this case is 
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Figure 16. Time- averaged kinetic energy densities and Nusselt number Nm of non- magnetic 
convection (thick lines) and of quadrupolar (diamonds), mixed (stars), hemispherical (triangles) 
and chaotic dipolar dynamos (squares). Both non-magnetic convection and dynamo solutions 
are in the case P = 1 and r = _10 4 . Values of P m of the dynamo cases decrease from 10 to 0.4 
as R increases. The densities Et, Et and Nth are shown in the left, middle and right panel, 
respectively. 



apparent from figure 15. The axisymmetric components of the magnetic field are almost 
steady in time. But the differential rotation in the polar regions may change its sign. 
At the particular time of figure 15 the zonal flows near the north- and south-poles have 
opposite signs. The fluctuating component of convection is dominated by the m = 1-mode 
and exhibits the attachment to the equatorial boundary as must be expected for inertial 
convection according to the discussion in section 3. For Prandtl numbers of the order 5 
and higher, on the other hand, the energy of the mean poloidal field always exceeds that 
of the mean toroidal field, since the differential rotation which creates the latter is not 
strong enough. But this changeover is rather gradual and depends more on the Rayleigh 
number than the other changeover discussed above. 

Since in some of the cases of figure HH the kinetic energies of convection in the absence 
of a magnetic field are shown, the effect of the magnetic field on convection becomes 
evident. In particular in the cases with P of the order unity or less the energy of the 
mean toroidal velocity field is strongly reduced by the dynamo action in that the Lorentz 
force leads to a braking of the differential rotation as has been discussed by Grote and 
Busse (2001). On the other hand the fluctuating components of convection are enhanced 
since the magnetic field prevents the onset of relaxation oscillations in which convection 
occurs only intermittently as has been discussed in connection with figure [2J This effect 
is clearly demonstrated in case (b) of figure [14] where the results for convection with 
and without magnetic field are compared for R — 2 x 10 6 . The same features are also 
revealed in figure [TH] where a more general comparison between convection with and 
without magnetic fields is shown. This figure also demonstrates that the symmetry of 
the magnetic field does not play a major role in the considerations of energetic aspects 
of dynamos. 

For Prandtl numbers larger than unity the effect of the magnetic field is much reduced 
which is evident from the fact that the increasing magnetic energy for increasing P m 
does not change the energy of the fluctuating components of convection as can be seen 
in figure [13] Only the differential rotation is reduced through the strengthening of the 
Lorentz force. The diminished influence of magnetic field on convection is also evident 
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from the property that typically viscous dissipation exceeds Ohmic dissipation by far for 
values of P of the order 5 or higher, while the two sinks of energy are more comparable 
at low values of P. 

6. Validity of the magnetostrophic approximation 

Because convection driven dynamos in rotating systems depend on a rather large num- 
ber of parameters it is desirable to eliminate one or more parameters through reductions 
of the basic equations. Among the nonlinear advection terms the momentum advection 
term appears to be most expendable since it does not seem to be essential for convection 
driven dynamos. In this way the magnetostrophic approximation is obtained in which 
the acceleration of fluid particles is neglected in comparison to the Coriolis force and 
the Lorentz force. In general this approximation can be easily obtained through the re- 
placement of the viscous time scale e? 2 jv by the time scale d 2 / (k 1 ^ 1 X 1 ') with ^ 7 ^ 1 
which is intermediate between the thermal and magnetic diffusion time scales. When 
y/ p[iK Y ~~t X< v I d is used as scale of the magnetic field the basic dimensionless equations 
of motion, equation of induction and the heat equation can be written in the form 

(pi-7p7)-i(a tM + u . v«) + T k x u = -Vtt + 9r + V 2 m + (VxB)xB (6.1a) 

{Pm/Py (dtB + u ■ V.B — B ■ Vu) = V 2 B (6.16) 

{p/Pra) 1 ^ (<9 t e + u ■ ve) = Ru r + v 2 e (e.ic) 

where k is the unit vector parallel to the axis of rotation. From the form of equations 
(|6.1jl it is clear that the magnetostrophic approximation should certainly be valid in the 
limit P 1 ~ 1 P^ l — > 00. In the following we shall focus on the case 7 = 0. In figure IT7l 
the energy densities have been plotted for fixed values of R,t and of ft/A. Since the 
fluctuating poloidal energy density E p always amounts to about 50 percent of the cor- 
responding toroidal one it has not been plotted. It can be seen that the kinetic energies 
tend to become independent of P with increasing P in accordance with the magne- 
tostrophic assumption. The energy density E t representing the differential rotation is 
the only exception as expected. Little indication of an approach towards the validity of 
the magnetostrophic approximation is found, however, when the magnetic energy densi- 
ties are considered. As has already been mentioned the dynamo process is rather sensitive 
to the presence of the differential rotation and much higher values of P may be needed 
before the magnetostrophic regime is approached. It is remarkable to see the distinct 
minimum of magnetic energies near P = 8 which corresponds to change in the structure 
of the magnetic field as has been mentioned. The transition from a geostrophic differential 
rotation to a thermal wind one appears to be mainly responsible for this feature. 

Results obtained on the basis of the magnetostrophic approximation in the case 7 = 
are by definition independent of P. In particular the ratio between the magnetic energy 
and kinetic energy will be proportional to P (Glatzmaier and Roberts, 1995) as is evident 
from the different scales used for the velocity and for B/^fpji in equation (|6.1a|) . The 
ratio between Ohmic and viscous dissipation would be independent of P and would 
depend only on n/\. This latter parameter seems to be even more important than the 
magnetic Prandtl number for convection driven dynamos. At least for Prandtl number 
of the order unity or less the condition 1 < P m /P appears to apply for the minimum 
value P m for which dynamo action can be obtained as is indicated by a comparison of 
dynamos obtained for P = 1 (see figure 1 of Grote et at, 2002) and for P = 0.1 (see 
figure [9]) . 

Similar conclusions as have been reached for 7 = may be reached for 7 = 1. In this 
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Figure 17. Kinetic energy densities E p (in (a)), E t (in (bj) and _B t (in (c)), all multiplied 
by P 2 as functions of P in the case r = 5 x 10 3 . The dynamos corresponding to k/X = 1, 
R = 5 x 10 5 are indicated by solid squares, to k/X = 1, R = 6 x 10 5 by solid triangles, to 
k/A = 1, it = 8 X 10 5 by solid circles, to k/X = 2, R = 6 x 10 s by crosses, to k/A = 2, _R = 10 6 
by plus-signs, to k/X — 5, R = 6 x 10 5 by empty squares, to k/A = 0.5, R = 10 6 by empty 
triangles and to non-magnetic convection with R = 10 6 by empty circles. The second row shows 
the corresponding magnetic energy densities multiplied by P. 



case the energy densities of figure [T7] can be plotted as a function of P m . For the cases 
k/X = 1 the plots remain unchanged with P m replacing P at the abscissa. For other 
values of k/X the respective data would just be shifted. Since no particular new insight 
appears to be evident we have not included such a figure. The general conclusion to be 
drawn is that the magnetostrophic approximation seems to become valid only for fairly 
high values of the parameter P 1_7 P^ with ^ 7 ^ 1. It should also be remembered 
that properties such as inertial convection and geostrophic differential rotation generated 
by Reynolds stresses disappear when the magnetostrophic approximation is used. Kuang 
and Bloxham (1996) have introduced for this reason a compromise in that the magne- 
tostrophic approximation is applied only to the non-axisymmetric part of the equations 
of motion. But their main result that the ratio of magnetic to kinetic energy reaches 10 3 
for P = 1 is far removed from properties of other convection driven spherical dynamos 
described in the literature. 



7. Effects of boundary conditions on dynamo solutions 

Because published solutions of convection driven dynamos in rotating spherical shells 
are based on a variety of parameter values and satisfy different mechanical, thermal 
and magnetic boundary conditions, a systematic comparison between various models is 
difficult to conduct. An attempt to investigate the differences between numerical models 
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Table 1. Time-averaged global properties of dynamos with various velocity and magnetic 
boundary conditions as follows. A: stress-free and insulating, B: no-slip and insulating, C: 
no-slip and a finitely-conducting inner core and D: stress- free and a perfectly-conducting inner 
core. The predominant symmetry type is indicated with "D" if dipolar, "Q" if quadrupolar and 
"— " if the dynamo is decaying. In the case B in left part of the table results for both a dynamo 
solution and a non-magnetic convection case are given. 
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based on various assumptions has most recently been made by Kutzner and Christensen 
(2000, 2002) who consider thermal and compositional driving with internal or external 
distribution of energy sources and a variety of thermal boundary conditions. Here we 
complement their studies with a few remarks on the possible effects of various velocity 
and magnetic boundary conditions. For this purpose we have selected two typical dynamo 
solutions. The dynamo at parameter values P — 0.1, r = 10 5 , R = 4 x 10 6 , P m = 0.5 
has a dipolar symmetry while the solution at P — 5, r = 5 x 10 3 , i? = 8x 10 5 , P m = 3 
has a quadrupolar one. The results for the standard boundary conditions (|2.7p and (|2.8p 
used in this paper are listed under A. We have repeated the simulations of these two 
cases for three of models with different combinations of boundary conditions and present 
the results in table [TJ For type B no-slip and insulating conditions have been adopted, 
for type C we have chosen no-slip conditions and a finitely-conducting inner core and 
finally type D has stress-free velocity boundary conditions and a perfectly-conducting 
inner core. Thermal conditions have remained unchanged in that temperatures are fixed 
at the boundaries in all four models. Snapshots of the spatial structure of the magnetic 
field are shown in figures IT51 and [TTJl 

The replacement of the stress-free condition by the no-slip condition at the inner and 
outer boundaries leads to considerable changes in the structure of the velocity held which 
are even more dramatic in the absence of a magnetic field. The differential rotation is 
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Figure 18. Effects of various boundary conditions on dynamo solutions for P = 0.1, r = 10 5 , 
R = 4 x 10 6 , P m = 0.5 in the dipolar cases A, B, C (from left to right) of table [T] The left 
halves of each plot show lines of constant zonal magnetic field B v and the right halves show 
meridional field lines, rsmOdeh = const., at particular moments in time. 




Figure 19. Same as figure [18] but for the quadrupolar cases D and A and D for the parameter 
values P = 5, r = 5 x 10 3 , R = 8x 10 5 , P m = 3. 

strongly inhibited by the Ekman layers which tend to oppose any deviation from rigid 
body rotation. The circulation induced by the Ekman layers, on the other hand, promotes 
the heat transport and and is responsible for the fact that the energy component E p is 
typically much larger in the case of no-slip boundaries than in the case of stress-free 
ones. But the main effect of the rigid boundaries is the suppression of the destruction 
of convection columns by the shearing action of the differential rotation. Relaxation 
oscillations which are predominant for Prandtl numbers of the order unity or less in 
the presence of stress-free conditions have not yet been found when no-slip conditions 
are used. It must be expected, however, that the properties of convection with rigid 
boundaries will approach those of convection with free boundaries as the Coriolis number 
t increases and the influence of the Ekman layer decreases. 

The effects caused by the differences in the velocity boundary conditions are reduced 
in the presence of a magnetic field generated by the dynamo process. The differential 
rotation is reduced by the action of the Lorentz force for both types of boundary con- 
ditions. But while the energy of the fluctuating components of convection is amplified 
by the presence of the magnetic field this effect usually does not happen in the case of 
no-slip boundaries as is evident from case B of table [TJ On the contrary, the magnetic 
degrees of freedom compete with kinetic ones for the same source of energy. Nevertheless, 
in spite of the lower kinetic energies in case B in comparison to those of case A, case B 
still exhibits a higher convective heat transport as is evident from the Nusselt number. 

The comparison of convection driven dynamos in case of P = 5 as presented by the 
right side of table [T] has been less successful since it was not possible to obtain dynamos 
in the cases B and D even after the Rayleigh number has been increased up to twice 
the indicated value. Except for the axisymmetric components of the velocities the kinetic 
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energies of the cases A and B are surprisingly close which indicates that the effect of 
no-slip boundaries on convection is similar to that of a dynamo generated magnetic field. 
The closeness of the energy results of the cases A and D for both values of P and of cases 
B and C for P = 0.1 indicates that the magnetic boundary conditions at the inner core 
have rather little influence on the dynamo process. It is remarkable, however, that in case 
D a quadrupolar dynamo is realized instead of the dipolar one in the standard case A. 
Even more remarkable is the property that this quadrupolar dynamo does not oscillate in 
time as has been found for all quadrupolar dynamos with the boundary conditions (|2.8p . 
From a general point of view it must be expected that for turbulent dynamos the inner 
core can not exert a strong influence on the dynamo process since the highly fluctuating 
magnetic field will not penetrate far into the inner core nor is it likely that strong fluxes 
will accumulate there since they are subject to Ohmic decay. Wicht (2002) has arrived 
at similar conclusions. 



8. Concluding remarks 

In this paper the analysis of convection driven dynamos has focused on the Prandtl 
number dependence while other parameters such as R and r have not been varied as 
much, but instead have been kept at reasonably high values which are numerically acces- 
sible with adequate resolution in space and in time. In the case of the Rayleigh number 
the numerical limitations are not felt as strongly as in the case of r since dynamos often 
seem to disappear or at least exhibit decreasing energies of the axisymmctric components 
of their magnetic fields as R is increased beyond an optimal value. More restrictive are 
the limits on the Coriolis number r which have prevented so far the attainment of dy- 
namos dominated by Ohmic dissipation. While the latter exceeds the viscous dissipation 
in some of the low Prandtl number dynamos, sufficient viscous friction always seems to 
be required to obtain numerically well-behaved solutions. Ohmic dissipation also helps 
in this respect and permits solutions at values of R and r which are not accessible at 
the same numerical resolution without magnetic field. But Ohmic dissipation does not 
necessarily prevent tangential discontinuities of the velocity field which are possible in 
the absence of viscous friction. 

A primary goal of dynamo simulations are results which can be compared quantita- 
tively with properties of the geomagnetic field and its variations in time. The geodynamo 
appears to exhibit features which are similar to those of convection driven dynamos with 
large values of P as well as to those with P of the order of unity. The fact that the mag- 
netic energy exceeds the kinetic energy in the Earth's core by a factor of the order 10 3 
together with the property that the axisymmetric poloidal component of the geomagnetic 
field dominates over the non-axisymmetric components suggest that the geodynamo re- 
sembles a high-P dynamo. On the other hand, the strong variations of the amplitude of 
the magnetic field on the magnetic diffusion timescale together with the torsional oscil- 
lations which manifest themselves as "jerks" on a much shorter timescale (Bloxham et 
ah, 2002) indicate a relationship to dynamos with P < 1 (Busse and Simitev, 2004b). Of 
course, these similarities may change when more realistic dynamo simulations with much 
lower values of the magnetic Prandtl number P m become possible. But they could also 
indicate that the two sources of buoyancy apparently present in the Earth's core, - the 
thermal one characterized by P < 1 and the compositional one characterized by P 3> 1 , 
should be distinguished in dynamo simulations. An indication that new dynamical effects 
become important in the presence of buoyancy sources with rather different diffusivities 
has been given by the linear analysis of Busse (2002b). 
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